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We discuss how the presence of gauge subsystems in the Bacon-Shor code [D. Bacon, Phys. 
Rev. A 73, 012340 (2006)] leads to remarkably simple and efficient methods for fault-tolerant error 
correction (FTEC). Most notably, FTEC does not require entangled ancillary states and it can be 
implemented with nearest-neighbor two-qubit measurements. By using these methods, we prove a 
lower bound on the quantum accuracy threshold, 1.94xl0 -4 for adversarial stochastic noise, that 
improves previous lower bounds by nearly an order of magnitude. 
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Operating a full-scale quantum computer will require 
methods for protecting against decoherence or system- 
atic hardware imperfections. One of the central results 
in the theory of fault-tolerant quantum computation [l[ , 
the quantum threshold theorem shows that a noisy quan- 
tum computer can accurately and efficiently simulate any 
ideal quantum computation provided that noise is local 
and its strength is below a critical value known as the 
quantum accuracy threshold 0]. The actual value of the 
accuracy threshold is therefore of great practical interest 
as it represents the desired target accuracy of prospec- 
tive implementations of quantum computation. However, 
determining this value is in general a very challenging 
problem since not only does it depend on the particular 
character of noise but it also depends on the specifics 
of the quantum computing architecture (e.g., geomet- 
ric constraints on qubit interactions). In this Letter, we 
present methods of fault-tolerant quantum computation 
based on a new subsystem code that significantly improve 
the existing lower bounds on the value of the quantum 
accuracy threshold and that are also beneficial for geo- 
metrically local quantum computing architectures. 

Subsystem codes protect quantum information from 
noise by mapping the system to be encoded (e.g., one 
logical qubit) into a subsystem of a larger system [H 0]; 
in fact, such a mapping provides the most general pos- 
sible encoding of quantum information More con- 
cretely, if we let TLs be the Hilbert space of the sys- 
tem to be encoded, the encoding maps density oper- 
ators and observables from Hs to a subsystem with 
Hilbert space Hl which lies inside a larger Hilbert space, 
H = {Ti-L® Ht) © Hr, where Ht describes additional 
"gauge" subsystems not necessarily protected from noise 
and Hr labels the "rest" of H. Fault-tolerant quan- 
tum computation has traditionally used subspace codes 
that can be seen as subsystem codes where Ht is one- 
dimensional encoding no subsystem. 

Interestingly, the converse is also true: For any sub- 
system code there exists a corresponding subspace code 
that can be obtained by "picking a gauge," i.e., by pro- 



jecting the state in Ht onto a pure state [4|]. Because 
of this connection, subsystem codes do not have proper- 
ties that make them fundamentally different from sub- 
space codes. Nevertheless, subsystem codes are in a cer- 
tain sense more efficient than the corresponding subspace 
codes since they require the extraction of fewer syndrome 
bits (but, see also [§]). As we will discuss in this Let- 
ter, subsystem codes are also advantageous for quantum 
computation because the presence of gauge subsystems 
can help us simplify the quantum circuits that implement 
fault-tolerant error correction. 

We will consider fault-tolerant quantum circuits where 
computation is encoded using a subsystem code due to 
Bacon [7| — because of the close relation of this code with 
Shor's code we will refer to it as the Bacon-Shor 
code. There is a different Bacon-Shor code for every in- 
teger n > 1; for fixed n, the corresponding code, C^J, 
is a distance-n stabilizer CSS code [l0( encoding one 
protected logical qubit into a code block of n 2 physical 
qubits. 

If we imagine placing the n 2 qubits in the C^g block 
on the vertices of an nxn square lattice as in Fig.[TJ the 
code's stabilizer group S, is generated by 
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where we have used the short-hand notation for the 
Pauli matrices A = cr x , Z = a Zl and Oj>, 0*j de- 
note operators that act nontrivially as a tensor product 
of O operators on all qubits in row or column j, respec- 
tively. Throughout, 1 m is understood to indicate the set 
{l,2,...,ro}. 

The code syndrome e — i.e., the binary vector of eigen- 
values of the 2(n— 1) stabilizer generators in Eq. (pj — 
induces a decomposition of the Hilbert space, H, of the 
n 2 qubits in the code block into subspaces encoding 
n 2 — 2{n— 1) = (n— l) 2 + 1 logical qubits. Therefore, 
within each subspace with fixed syndrome — and, in par- 
ticular, within the code space corresponding to the trivial 
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FIG. 1: Qubits in the C^J block sit on the vertices of an nxn 
square lattice. An element of the code stabilizer, the operator 
X2,*X3,* applies X on all qubits shown in blue. 



syndrome — we can define a subsystem decomposition 



(2) 



where we associate Hl with the one logical qubit pro- 
tected by the full distance, n, of the code. The logical 
Pauli operators for this logical qubit can be defined as 
Xl = X± .* (i.e., a tensor product of X operators applied 
on all qubits in the first row) and Zl = (i.e., a tensor 
product of Z operators applied on all qubits in the first 
column). The remaining (n— l) 2 logical qubits live in TCt 
and their logical Pauli operators can be chosen from the 
non-abelian group 

T = ( Xj^Xj+ij ; ZijZij + i \ieZ n ;j e 1 n -\ ) , (3) 

where we have used the notation Oij for an operator 
O acting on the qubit with coordinates Indeed, 
the operators in T commute with every operator in the 
code stabilizer, they commute with the logical operators 
Xl and Z^ and, furthermore, they can be grouped into 
(n— l) 2 independent pairs of anticommuting operators 
with operators in different pairs commuting. 

Given some nontrivial syndrome value, error recovery 
for the logical qubit encoded in rii, proceeds in a similar 
manner as in the classical repetition code: The eigenval- 
ues of the stabilizer generators {Xj^Xj+i.*} can be used 
to correct Z errors on up to rows. Moreover, only 

the parity of Z errors in each row is relevant: an operator 
acting as Z on a pair of qubits at the same row is an op- 
erator in T and, therefore, has no effect on the protected 
information in H.l ■ Error recovery for X errors proceeds 
similarly along the columns. 

Since the logical Pauli operators for the logical qubits 
encoded in TLt act nontrivially on only two qubits in the 
code block, if we take into consideration all (n— l) 2 + 1 
logical qubits then Cg'g has distance 2 and it is an error 
detecting code. However, if we only consider the logical 
qubit encoded in TLt then the effective distance is n and 
errors with support on up to L n /2J qubits in the code 



block can be corrected — we will call this logical qubit 
the protected qubit. In fact, error recovery for the pro- 
tected qubit may unavoidably result in applying at the 
same time nontrivial logical operations with support on 
TLt- This is not a problem as long as we never encode 
quantum information in TLt- We can think of the (n— l) 2 
logical qubits encoded in TLt as gauge qubits since they 
correspond to degrees of freedom for the logical informa- 
tion encoded in TLl- In some cases it will be sufficient to 
completely disregard the state of the gauge qubits. More 
interestingly, we will next discuss how fault-tolerant er- 
ror correction for the protected qubit can be simplified 
by taking advantage of the presence of gauge qubits. 

Let us first explain how we can extract the code syn- 
drome indirectly by manipulating the state of the gauge 
qubits. Our first observation is that we can express the 
stabilizer generators in Eq. (TT]) as 
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What is remarkable about this decomposition is that the 
operators in parentheses are supported on TLt] hence, 
they commute with operators in the code stabilizer and 
also commute with the logical operators for the protected 
qubit. Because of this, we can measure each of them sep- 
arately and then Eq. ^ implies that the code syndrome 
can be computed by taking the appropriate parities of the 
measurement outcomes. Moreover, since these operators 
act nontrivially on only two qubits in the code block, 
measuring them is especially easy; e.g., Fig. [2] shows sim- 
ple circuits for measuring Xj^Xj+i^ and ZkjZ, 
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FIG. 2: (a) A circuit for measuring the operator Xj^Xj+i^ 
using one ancillary qubit. (b) A similar circuit for measuring 
Zt,jZk,j+i- |+) oc |0) + |1) is the +1 eigenstate of X. 

This indirect method for inferring the code syndrome 
significantly reduces the overhead for fault-tolerant error 
correction (FTEC). This is because, unlike in standard 
FTEC methods [|J , this method does not require prepar- 
ing and verifying entangled ancillary states — ancillary 
qubits in the |0) or |+) state are sufficient as in Fig. [2] 
Furthermore, for maximum qubit efficiency but at a cost 
of extra memory error, even a single ancillary qubit suf- 
fices to sequentially measure all two-qubit operators nec- 
essary to extract the syndrome! In addition, the specific 
two-qubit operators to be measured can be chosen to have 
support on neighboring qubits in the code block when 
these qubits are arranged on a two-dimensional square 
lattice. For this reason, this FTEC method could prove 
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to be especially advantageous for quantum computing 
with geometric locality constraints. 

In settings where nonlocal qubit interactions are pos- 
sible or when qubit movement has error rate much lower 
than quantum gates, other standard FTEC methods will 
yield the best accuracy thresholds. Since C^g is a CSS 
code, Steane's FTEC method [l2| can be used to extract 
the syndrome provided we can fault-tolerantly prepare 
logical |0) (|0)i) and logical |+) states for the 

protected qubit. Alternatively, KnilPs FTEC method 
[131 ] requires fault-tolerantly preparing two code blocks 
where the two protected qubits are encoded in a logical 
Bell state (|$ }l oc |0} l |0) l + 

As is evident from the decomposition in Eq. ([2]), the 
distinctive feature of subsystem codes is that logical 
states in Tij, are not uniquely encoded: after having spec- 
ified a logical state in Hl, the state in Ttr can still be 
arbitrary. This freedom in choosing the state of the gauge 
qubits can be used to our advantage in the design of en- 
coding circuits of logical states for the protected qubits. 
In particular, we will next discuss how we can design re- 
markably simple encoding circuits for the logical ancillary 
states required for Steane's and Knill's FTEC methods, 
thus also reducing the overhead associated with the post- 
encoding verification of these states. 

For concreteness, consider designing an encoding cir- 
cuit for |0}l, i.e., a state in the code space which is the 
+1 eigenstate of Zl- With the state in TLl specified, we 
can choose the state in Ht to be the +1 eigenstate of 
the gauge-qubit operators {Xi,jXi + i,j}. In other words, 
our encoding circuit prepares the +1 eigenstate of the 
operators in the following stabilizer group: 
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We recognize the state described by Eq. ([5]) as a tensor 
product of n Schrodinger "cat" states in the Hadamard- 
rotated basis, each one lying across a column in Fig.Q] 

We can obtain the state by applying a logical 

Hadamard transformation to the state \0)l- We observe 
that applying Hadamard gates bitwise has the same effect 
as a logical Hadamard gate up to a 90-degree rotation of 
the square lattice in Fig.[T] (and, also, up to a nontrivial 
operation acting on the gauge qubits). We can therefore 
obtain by first preparing \0}l, applying Hadamard 
gates bitwise and, finally, rotating the lattice in Fig. Q] by 
90 degrees. The bitwise Hadamard gates will produce 
usual cat states in the computation basis which, after 
rotating the lattice by 90 degrees, will be aligned each 
to lie across a row — this is our construction for \+)l- 
In addition, logical Bell states for Knill's FTEC method 
can be constructed by interacting two blocks encoded in 
the states \+)l and |0}l via a logical CNOT gate. Since 
the logical CNOT gate can be implemented by transversal 
CNOT gates, this construction is also especially simple. 

The logical CNOT and Hadamard gates have simple 
transversal implementations and, together with single- 



qubit preparation and measurement, they suffice for per- 
forming fault-tolerant error correction; however, they are 
not universal for quantum computation. Universal quan- 
tum computation can be realized by including the logical 
phase gate, S = exp(— ijcr z ), and one logical non-Clifford 
gate in our gate set; for a detailed discussion on achieving 
encoded quantum universality see chapter 5 in [111 ]. 

We have analysed the accuracy threshold of fault- 
tolerant quantum circuit simulations that use the con- 
catenated Cgg for various values of n. In our analysis, 
we have considered adversarial stochastic noise 
the stochastic form of local noise [111. ll4Lll5l] . In this noise 
model, we assume that any r specific elementary physi- 
cal operations in the noisy quantum circuit (single-qubit 
preparations, quantum gates, memory steps or single- 
qubit measurements) fail with probability at most e r for 
some constant < e < 1. Noise is adversarial because 
faults need not be independent and, moreover, noise may 
act jointly on the support of all different faulty opera- 
tions. Adversarial stochastic noise is the natural form 
of noise to consider when analyzing fault-tolerant circuit 
simulations that use concatenated codes. This is because 
even if physical noise does not include adversarial corre- 
lations between different faulty locations, the effective 
noise that acts on all coding levels of the concatenation 
hierarchy higher than the physical level may include such 
correlations which arise due to coding and the propaga- 
tion of syndrome information [Til ] . 

Our rigorous lower bounds on the quantum accuracy 
threshold were obtained by performing an analysis of ma- 
lignant sets of locations on extended rectangles accord- 
ing to the method that was introduced in [14[. Table 
U summarizes our results. We carried out the required 

(3) 

combinatorial analysis for the concatenated C BS using 
both Steane's and Knill's FTEC methods and for the 

(5) 

concatenated C BS using only the former method due to 
time limitations. The analysis was done by using a com- 
puter program running for several months on a cluster 
of 20 Pentium III processors. Our best rigorous lower 
bound on the accuracy threshold, 1.94 x 10~ 4 , was ob- 

(5) 

tained for the concatenated C BS . This lower bound im- 
proves by nearly an order of magnitude the 2.73 xlO -5 
lower bound that was proved in [14J with the concate- 
nated Steane [[7,1,3]] code. 

(5) 

Analyzing codes of larger block size than C BS proved to 
be computationally difficult in this exact setting. In these 
cases, we have used a Monte-Carlo method to uniformly 
sample the set of fault paths with a fixed number of faulty 
locations inside an extended rectangle. By estimating 
what fraction, /, of these sets is malignant, we obtain 
an estimate of the exact combinatorial coefficients that 
determine the accuracy threshold with a standard error 



that can be determined by using a j = y f(l — f)/N 
where N is the sample size. We have also applied this 
Monte-Carlo method to the cases where we could ex- 
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Steane 
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Knill 


2,681 




1.91 ±0.01 



TABLE I: Rigorous lower bounds on the accuracy thresh- 
old, eo, for adversarial stochastic noise with the concatenated 
Bacon-Shor code of varying block size and comparison with 
prior rigorous lower bounds using the concatenated Steane 
[[7,1,3]] code and Golay [[23,1,7]] code 0. The third 
column gives the number of locations in the CNOT extended 
rectangle [l4|. The forth column gives exact lower bounds on 
£o; the results are obtained using a computer-assisted com- 
binatorial analysis. The fifth column is the Monte-Carlo es- 
timate for eo with la uncertainties. Bold fonts indicate the 
best results in each column. 



tract the exact combinatorics in order to provide evidence 
that the Monte-Carlo estimates are accurate. Indeed, as 
can be seen in Table U the exact lower bounds in those 
cases lie within la of the estimated lower bounds. Our 
Monte-Carlo results give evidence that the quantum ac- 
curacy threshold achieves a maximum over all concate- 
nated Bacon-Shor codes of varying block size for the 25- 
qubit code using Knill's FTEC method. Chapter 5 in 
[111 ] discusses the details of our threshold analysis for the 
CNOT extended rectangle and explains why the accuracy 
threshold is determined by this gate and not by non- 
Clifford gates. 

In conclusion, we have shown how the presence of 
gauge qubits in the Bacon-Shor code can be exploited 
to design quantum circuits for fault-tolerant error cor- 
rection with remarkable properties. We have presented a 
new method for fault-tolerant error correction that uses 
nearest-neighbor two-qubit measurements and does not 
require the preparation of entangled ancillary states. We 
expect this method to be advantageous in implementa- 
tions that impose geometric locality constraints in the 
interactions between qubits such as, e.g., in proposed 
ion-trap or solid-state schemes of quantum computa- 
tion. On the other hand, standard methods for fault- 
tolerant error correction can be implemented using espe- 
cially simple encoding circuits for the required ancillary 
logical states, thus greatly reducing the overhead asso- 
ciated with the verification of these states. Our lower 
bound on the quantum accuracy threshold, 1.94x10 
for adversarial stochastic noise, is the best that has been 
rigorously proven so far and improves previous rigorous 



lower bounds by nearly an order of magnitude. 

An open challenge for future work is to investigate to 
what extent the lower bounds we have presented can be 
improved by performing a more detailed combinatorial 
analysis or by considering a less adversarial noise model. 
In addition, it would be important to derive lower bounds 
on the accuracy threshold in a setting where qubit inter- 
actions are geometrically local in order to quantify the 
advantages of our FTEC methods and to compare with 
other fault-tolerant schemes for local quantum computa- 
tion such as, e.g., those in 17, 18j |. 
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